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2 ■ Abstract 
O 

. We study the dynamical evolution of a stationary, axisymmetric, and perfectly conducting cold 



accretion disk containing a large-scale magnetic field around a Kerr black hole, trying to under- 
stand the relation between accretion and the transportation of angular momentum and energy. A 
■ one-dimensional radial momentum equation is derived near the equatorial plane, which has one 

, intrinsic singularity at the fast critical point. We solve the radial momentum equation for solutions 

> 

Cf^ ' corresponding to an accretion flow that starts from a subsonic state at infinity, smoothly passes the 

o ■ 

' fast critical point, then supersonically falls into the horizon of the black hole. The solutions always 

(N 



O 



have the following features: 1) The specific energy of fluid particles remains constant but the spe- 
cific angular momentum is effectively removed by the magnetic field. 2) At large radii, where the 



Qh' disk motion is dominantly rotational, the energy density of the magnetic field is equipartitioned 



with the rotational energy density of the disk. 3) Inside the fast critical point, where radial motion 



o 

\ becomes important, the ratio of the electromagnetic energy density to the kinetic energy density 

drops quickly. The results indicate that: 1) Disk accretion does not necessarily imply energy dissi- 

X 

■ pation since magnetic fields do not have to transport or dissipate a lot of energy as they effectively 

transport angular momentum. 2) When resistivity is small, the large-scale magnetic field is ampli- 
fied by the shearing rotation of the disk until the magnetic energy density is equipartitioned with 
the rotational energy density, ending up with a geometrically thick disk. This is in contrast with 
the evolution of small-scale magnetic fields where if the resistivity is nonzero the magnetic energy 
density is likely to be equipartitioned with the kinetic energy density associated with local random 
motions (e.g., turbulence), making a thin Keplerian disk possible. 

PACS numbers: 04.70.-s, 97.60.Lf 
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I. INTRODUCTION 



Accretion onto a gravitating object is believed to be a principal energy source for powering 
many astrophysical systems, including active galactic nuclei (AGN), quasars, and some 
types of stellar binaries [1]. Since accreting gases far from the central object typically have 
sufficiently large specific angular momentum compared to a particle moving on a circular 
orbit around and near the central object, the formation of an accretion disk is unavoidable. 
To accrete onto the central object and release gravitational energy, gases at large distance 
must lose a lot of angular momentum. It is usually assumed that some kind of viscosity 
slowly operates in a disk to transport angular momentum outward and dissipate energy, so 
that to a very good approximation in the disk region particles move on circular orbits with 
a small radial velocity superposed. 

The simplest stationary and axisymmetric model of accretion disk is a geometrically thin 
Keplerian disk, where it is assumed that in the disk region the gravity of the central object 
is approximately balanced by the centrifugal force so that disk particles move approximately 
on Keplerian circular orbits 0,0,0]. Then, at each radius, the specific angular momentum 
and the specific energy of disk particles are given by that of a test particle moving on the 
Keplerian circular orbit at that radius. When the central object is a black hole, the Keplerian 
disk has an inner boundary at the marginally stable circular orbit, inside which the gravity is 
too strong to be balanced by the centrifugal force 0, 0, . The total efficiency of converting 
mass into energy is defined by 

6S 



6M 



(1) 



where 6S is the total released energy during the process of accretion, 6M is the corresponding 
accreted mass 0|. Then, the total efficiency of a Keplerian disk around a Kerr black hole 
is given by the difference between the specific energy at infinity {Eoo = 1) and the specific 
energy on the marginally stable circular orbit (-Ems) 

^Keplerian 1 ^ms • (2) 

Thin Keplerian disks correspond to the most efficient case of releasing the gravitational 
energy of an accretion fiow. For a thin Keplerian disk around a Schwarzschild black hole, 
the total efficiency is ^Keplerian ~ 0.06. For a thin Keplerian disk around a "canonical" Kerr 
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black hole of specific spinning angular momentum a = 0.998M, where M is the mass of 
the black hole, the total efficiency is exepierian ~ 0.32 when the effect of radiation capture 
is considered p]. These efficiencies of converting mass into energy are indeed much higher 
than nuclear fusion, the latter is ~ 0.007 for hydrogen. The released gravitational energy 
in a thin Keplerian disk is assumed to be efficiently dissipated by the internal viscosity and 
instantly radiated away from the surfaces of the disk. 

Though the model of geometrically thin and optically thick Keplerian accretion disks is 
supported by the appearance of Big Blue Bump in the spectra of all bright AGN 
and references therein), it is challenged by the observational fact that most nearby galactic 
nuclei are much less active — sometimes not active at all and references therein). These 
objects have luminosities that are lower than what a thin Keplerian disk model usually 
predicts by many orders of magnitudes. In order to interpret this phenomenon, different 
disk models are proposed where either the radiative efficiency is assumed to be low so that 
the released gravitational energy is converted into entropy and internal energy trapped in 
the flow and carried into the central black hole by the flow, or the mass accretion rate is 
assumed to be low and most of the released gravitational energy is carried away by a mass 
outflow 0, 12, 13j_14, 15, These alternative models are summarized by Blandford 
jl7 | and Naravan |lol|. Though they predict a different disk luminosity, these models are 
in common with the thin Keplerian disk in the following aspect: the mechanical energy (= 
gravitational potential energy + kinetic energy) of the disk flow is efficiently converted into 
other forms of energy. 

In the theory aspect, the proposed disk viscosity is poorly understood. Investigations on 
the magnetohydrodynamics (MHD) in a differentially rotating disk suggest that magnetic 
fields may play the role of disk viscosity ( 3, and references therein) . Due to the exis- 
tence of magnetorotational instability, MHD turbulence is easily generated in a differentially 



rotating dis 



which effectively transports angular momentum outward and drives accretion 



, very few models of disk radiation 
are based on these ideas, and generally angular momentum transport is treated through a 
parameterized viscous stress — the so-called a- viscosity following Shakura & Sunyaev 2, Ol ■ 
This is because of the fact that though it is well established that magnetic fields effectively 
transport angular momentum and drive accretion, it is far from clear how magnetic fields 
transport and dissipate energy in the meantime. Numerical simulations are a powerful tool 
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for studying the nonlinear evolution of magnetic fields. However, numerical simulations usu- 
ally have large numerical dissipation due to finite spatial resolutions, which makes energy 
transportation and dissipation cannot be accurately studied with current numerical simula- 



tions |l8, 



IS 



25 



26|. Therefore, our current understanding on how magnetic fields transport 



and dissipate energy is much worse than our understanding on how magnetic fields transport 
angular momentum. 

In this paper we ask and try to answer the following question: Do magnetic fields transport 
or dissipate energy as efficiently as they transport angular momentum? Or, more precisely, 
Do magnetic fields have to transport or dissipate a lot of energy as they efficiently transport 
angular momentum in a disk accretion flow? 

To answer this question (or, to get some insight into the answer), we use an analytical 
model to study the dynamical effects of magnetic fields in an accretion disk around a Kerr 
black hole. The disk is assumed to be stationary, axisymmetric, perfectly conducting, and 
contain a large-scale magnetic field j^. Both the magnetic field and the velocity field 
of the disk fluid have only radial and azimuthal components in a small neighborhood of 
the disk central plane 9 = 7r/2, where 6 is the polar angle from the symmetry axis of the 
black hole. In other words, near the central plane of the disk, each magnetic field line 
lies on a surface of constant 6, and the same is true for each fluid stream line. (Similar 
approach is usually used to study the infiow/outfiow problem associated with a low densit\ 



23, m. 



)ciated witn a low den sity 
and 13311 



magnetosphere around a black hole or neutron star, see 
for a review.) For simplicity, we neglect the thermal pressure in the disk by assuming that 
the accretion flow is cold. Then, making use of the conservation of rest mass, energy, and 
angular momentum. Maxwell's equations and dynamical equations are reduced to a set of 
one-dimensional equations with radius r as the only variable. We will self-consistently solve 
these equations, to look for the relation between accretion and the transportation of angular 
momentum and energy. 

Gammie Q] has used the model to study the flow on the equatorial plane in the plunging 
region between the inner boundary of the disk and the horizon of the black hole. Gammie's 
model is further explored by Li |35| and tested with numerical simulations by De Villiers 
& Hawley 3^. In this paper we extend the model into the disk region, and into a small 
neighborhood of the equatorial plane. This extension allows us to investigate a disk accretion 
flow driven by magnetic flelds from a global view: from the subsonic far outer disk region 
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down to the supersonic deep plunging region 



3- 



We emphasize that such a model can only 



be applied to a large-scale and ordered magnetic field. For small-scale and chaotic magnetic 
fields, the gradient of the magnetic field is important, making the solutions invalid when 
even a small nonzero resistivity is introduced (see Sec. I VII for detail). 

The paper is organized as follows: In Sec. [H] we exactly solve Maxwell's equations for 
the model outlined above. In Sec. II I II we write down the equations for the conservation of 
rest mass, the conservation of angular momentum, and the conservation of energy. We show 
that the solutions always correspond to an accretion fiow with constant specific energy when 
reasonable boundary conditions at infinity are satisfied. Then we derive a one-dimensional 
radial momentum equation corresponding to the magnetized disk accretion fiow. In Sec. IIVI 
we discuss the analytical properties of the radial momentum equation, the solutions on 
the horizon of the black hole, and the asymptotic solutions at infinity. We show that the 
fast critical point is the only intrinsic singularity in the radial momentum equation. In 
Sec. EJwe present some explicit solutions representing a magnetized disk fiow with constant 
specific energy accreting from infinity onto a central black hole. In Sec. IVII we discuss the 
effect of finite resistivity and compare the evolution of large-scale magnetic fields to that of 
small-scale magnetic fields. In Sec. I VI II we summarize the results and draw our conclusions. 

In Appendix^ we make some effort to generalize our results to more general and more 
complicated models. In Appendix El we present the Newtonian limit of our solutions. 

II. SOLUTIONS TO MAXWELL'S EQUATIONS 

We assume that the background spacetime is described by the metric of a Kerr black hole 

of mass M and specific angular momentum a, where —M < a < M . In Boyer-Lindquist 

coordinates (t, r, 9, 0), the metric of a Kerr black hole is given bv [btI. IssI ] 

, 2 / 2Mr\ ,2 4Mar . 2^ , ^ , 2 ^ ^ni Asin^^ , . 
ds^ = -il — dt^ — sin^ 9 dtdcj) + — dr^ + S d9^ + — - — d<p^ , (3) 

\ Z_( / Za 

where 

A = r^ -2Mr + a'^ , E = + cos^ 9 , A = {r^ + a^f ~ Aa^ sin^ 9 . (4) 

The radius of the outer event horizon is at rn = M + (M^ — a^)^^^. 

As stated in the Introduction, we focus on solutions in a small neighborhood of the 
equatorial plane defined by cos^ ^ 1. Then, to the first order of cos^ j57|, the metric does 



not depend on 6 and is given by 

(is^ = - 1 ]dt^ dtd4> + —dr^ + r^dO^ + — dc^\ 5 

where A = X{e = 7r/2). 

Assuming that the disk fluid is perfectly conducting so that magnetic field lines are frozen 
to the fluid. Then, in terms of coordinate components, the Maxwell equations that we need 
to solve are reduced to 

' ^-[^g{u-B^-u^B-)]=^^ (6) 



where x" = {t,r,6,(j)), \J —g = (to the first order of cos 6'), u" is the four- velocity of the 
fluid, and -B" is the magnetic field measured by an observer comoving with the fluid. (The 
electric field measured by an observer comoving with a perfectly conducting fluid is zero.) 
We assume that in the small neighborhood of the equatorial plane we have 

= , 5^ = , (7) 

i.e., each line of velocity or magnetic field lies on a surface of ^ = constant. Equation (jT)) 
implies that 

in the same neighborhood. Then, for stationary and axisymmetric solutions with d/dt = 
d/d(f) = 0, the Maxwell equation © is reduced to 

-^[r^u'B" ~u''B'')]=0. (9) 

By definition, B"^ satisfies UaB"" = 0. Then, making use of UaU"" = —1, we can solve for 
B^', B'f', and 5* from Eq. 0. The results are 

(10) 

(11) 
(12) 

where Cn and are constants. 



B' 




-CoUt + ^M^) , 


B^ 


^1 


[-CoMiM'^ + (l + MX)^] , 


B' 


1 
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The electromagnetic field tensor corresponding the solutions given by Eqs. ()10 p - ()12|) is 

Fab = eakcdu'B" = -2^ d^ade^ + ^ {C^u^ + ^u') dr^adOb] + 2Co de^adh] , (13) 
where eabcd is the totally antisymmetric tensor of the volume element that is associated with 



the metric 



38[, the brackets "[ ]" denote antisymmetrization. Gammie j3^ obtained the 
same solution for Fah by directly solving the Maxwell equation V[aFhc] = 0. 

It is useful to have a look at the magnetic field and the electric field measured by an 
observer in the locally nonrotating frame (LNRF observer, see P|) in the equatorial plane, 
whose four-velocity is f/" = [{9 / dtY + u {d / dcf))"'], where x = (r^A/A)^^^ and u = 
2Mar/A are respectively the lapse function (i.e., the redshift factor) and the frame dragging 
angular velocity. From the electromagnetic field tensor in Eq. (fT^ . we can calculate the 
magnetic field and the electric field measured by a LNRF observer 



2 bed ^2 ^ ^2^r 



{C,u^ + ^u^){£)\ (14) 



E'" - ^'^^^' = ^(^o- + v&)(|)' . (15) 

Therefore, the constant Co measures the flux of the radial magnetic field, and the constant 
\E' is related to the electric field measured by a LNRF observer and thus measures the 
direction of the magnetic field relative to the fluid motion js^. 

Equation (ITUI) can be rewritten as \E' = u'^^{r'^B^ + CoUt). For a Keplerian disk we have 
Ut = —1 and Utf, oc r^/^ as r ^ oo. If generally we assume that Ut and r'^B^ are bounded but 
Utf) is unbounded as r ^ oo, then we must have 

B'' 

^ = , i.e., — = — . (16) 

We note that it is not obvious that r'^B'^ must be finite, since the magnetic flux across a 
circle of a radius r is oc r^i?""/x, not r^i?*". Stronger justifications for the validity of Eq. ()16|) 
are given in the next section where dynamical equations are discussed. 
From Eqs. (fTn|l and (fTT|l we have 

B'i'u' ^ 
k ~ I = (17) 

B'u't' {-CoUt + ^M^) u't' ' ^ ^ 

Because of the appearance of u"^ in the denominator on the right-hand side, which generally 
quickly decreases with radius (e.g., u'^/u^ oc r^^/^ for a Keplerian disk), \k\ grows quickly 



as r increases, if \1/ 7^ 0. This means that, if at a large radius /c is ~ 1, then k quickly 
approaches zero as r decreases. 

As an example, assuming that at r = tq ^ rn we have k = k^ 1 (which is true if 
~ ri?*^ and \u^/ru'^\ ^ 1 at r = Tq), then by Eq. (jl7|) we have 

^ koUfU'^ 



f2vi/ = — 



r—ro 



(18) 



Co 1 + kou^u't' 

where f2o = {u^/u^^^^^ is the angular velocity of the disk at r = ro . In the last step we 
have used the fact that Ut ~ — m* ~ —1 and u^u'l' ~ M/r ^ 1 at r = tq . Substitute Eq. (fTH|) 
into Eq. (fTTj) . we then obtain 

{ut + koQoU^) u'f> 
at any radius r. For r -C ro in a Keplerian-type disk we have 



\k\ 



- « 1 , (19) 



since UqU^ ~ (M/ro)(r/ro)^^^ ^ 1 and ~ —1. 

When Eq. flTB|) is satisfied, from Eq. ([13)) we have E'"- = {Cquj/x't'^) {d/dOy. Then, for 
a Schwarzschild black hole (i.e., a = = 0), we have = 0, i.e., the magnetic field is 
parallel to the velocity field as seen by a local static observer (equivalent to a LNRF observer 
for a Kerr black hole). For a Kerr black hole with a 7^ 0, the electric field measured by a 
LNRF observer is nonzero. However, the electric field decays quickly with increasing r: 
E'^E'"" (X for r ^ th. Therefore, for a Kerr black hole, to a good approximation the 
magnetic field is also parallel to the velocity field in a LNRF, unless in the region close to 
the horizon of the black hole. Therefore, when Eq. (fTB|) is satisfied, we say that the magnetic 
field is parallel to the velocity field; or equivalently, the magnetic field lines are aligned with 
the fluid motion. But it must be kept in mind that this statement is not very accurate for 
the case of a Kerr black hole according to the above discussion. 

Equation (fT^ indicates that after the disk flow settles into the stationary and axisym- 
metric state the magnetic field and the velocity field tend to be aligned with each other if 
they are not so initially. 

In the next section we will see that a magnetic field satisfying Eq. (jlfjj) has a very in- 
teresting feature: it does not change the specific energy of fluid particles though it affects 
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the specific angular momentum. Physically this is caused by the fact that the electromag- 
netic force produced by the magnetic field satisfying Eq. |TBj) is always perpendicular to the 
velocity of the fluid. 



III. DERIVATION OF THE RADIAL MOMENTUM EQUATION 



For a cold and perfectly conducting magnetized fluid the total stress-energy tensor is 
given by T""^ = p^u'^u^ + T^m? where pm is the mass-energy density of the fluid matter as 
measured by an observer comoving with the fluid, is the stress-energy tensor of the 
magnetic field js^ 3^ 



An Svr An 

By Einstein's equation, the dynamical equations of the magnetized fluid are given by 



(20) 



. 



(21) 



Conservation laws 



35 



38, 



The con- 



A number of conservation laws can be derived from Eq. ()2ip 
traction of with Eq. (j^TJ) leads to the conservation of rest mass (also called the continuity 
equation): Va {pxaU"') = • The contraction of {d/d(f))°' with Eq. (j?T|) leads to the conserva- 
tion of angular momentum: VaT^" = . The contraction of {d/dtY with Eq. 1)2111 leads to 
the conservation of energy: VaT^"' = . [Note that {d/dtY and {d/d(f)Y ^^e Killing vectors 
of a Kerr spacetime.] 

Because of Eqs. ((Tj) and (jH)), for stationary and axisymmetric solutions in the small 
neighborhood of the disk central plane, the above conservation equations are simplified to 
one-dimensional ordinary differential equations js^ 

|: = ; (22) 

j-[r'[p^u^u' + T^^/)] = 0; (23) 
±[r^[p^,UtU^ + T^^Y)] = 0; (24) 

respectively. 



Substituting the solutions (fTIH) -(fT ^ into Eq. (pUjl . we obtain 

Tem/ = ^ {B'u.u^ - B,B^) = {Cou^ + ^u') , (25) 

^em/ = ^ [B'utu^ - B,B^) = -fi^T^M,/ , (26) 

where is defined by Eq. (ITH|) . (Though based on the arguments in Sec. HIl we expect 
\1/ = in a stationary state, in this section we keep ^ in the derived formulas to make the 
formulas as general as possible.) 
The solution to Eq. is 



= A ? r ' (27) 



where is the constant of radial mass fiux. From Eqs. pop - ()12|) . at large radii where 
r ^ M, Ut ~ — M* ~— 1, u'^^O, m''^~0, and is unbounded, we have = BaB"^ ^ 
\E'^/r^ (^'")^, if 7^ and u^/r^ <^ 1 (e.g., M^/r^ ^ M/r for a Keplerian disk). Then 

1 

^ oc , 28 

as r oo. Equation (j^Hj) implies that if \1/ 7^ 0, at sufficiently large radii the magnetic field 
must be so strong that B'^/Atc ^ p^. In order to prevent this to happen, \l/ must be zero. 
This provides a strong support for Eq. ()16|1 . 

Making use of Eqs. (j2SI)^(j2ZI), we can integrate Eqs. and to obtain 

+ ^ K - ^^^^0 = , (29) 
- (""^ - ^*"*) = , (30) 



where 



Co = ^^T- ' /L = > JE = , l-aij 



Fi^ = Attt'^TJ and F-^ = —Attt'^T^ ^ are the constants of radial angular momentum fiux and 
radial energy fiux, respectively. The dynamical equations and (jHn|l do not depend on 
the sign of cq. Without loss of generality, henceforth we assume that cq > . 
From Eqs. and (^01) we have 

E-n^L = -/e + , (32) 
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where E = —ut is the specific energy, L = -u,^ is the specific angular momentum of fiuid 
particles. Since /e, /l, and fi^ are constants, Eq. (jH^ immediately implies that the specific 
energy E must be constant when \1/ = (i.e., Vt^, = 0), and the specific angular momentum 
L must be constant when Cq = but ^ 7^ (i.e., = 00). 
The variation of Eq. (jH^ leads to 

6E = n^6L. (33) 

Obviously, Q^f, cannot be negative if we want 6E/6L > 0, i.e., angular momentum and 
energy propagate in the same direction. 

If the outer boundary of the disk fiow is at r = 00 and the specific angular momentum L 
is unbounded at r = 00, then we have 6L = L{r = 00) — L(r) = 00 for any finite r. Generally 
E is positive and bounded at r = cxd and E must be positive outside the ergosphere of a 
Kerr black hole. Hence, 6E = E{r = 00) — E{r) must be finite for r > 2M. Then, by 
Eq. (jnSl), this can be true if and only if Qq, = 0, which provides an additional support for 
Eq. ()16|) from the consideration of energy conservation. 

If the outer boundary of the disk fiow is at a finite radius tq ^ rn and Qi^, is given by 
Eq. (fTHj) . then, by Eq. (jH^ . for any r <^ tq we have 

6E = Eo- E{r) ^ k.U^L^ , (34) 

where Eq = E{r = vq), Lq = L{r = tq), and we have assumed that Lq ^ L{r). Notice that 
the value of 6E is determined by the boundary condition at the outer boundary of the disk 
when r <^ ro, which is in contrast with the case of a thin Keplerian disk where the total 
energy variation is determined by the boundary condition at the inner boundary of the disk. 
This is caused by the fact that in the present model the variation in energy can only happen 
in a region close to the outer boundary, since the magnetic field and the velocity of the fiuid 
quickly get aligned with each other in the region of r -C Tq. 

B. Dynamical equilibrium in the ^-direction 

Having solved the equations for the conservation of rest mass, angular momentum, and 
energy on each surface of ^ = constant near the equatorial plane, we need to check if the 
dynamical equilibrium in the ^-direction is maintained. 
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To do so, let us define a vector 

d_y ^ i_ 

dz ) sin V d9 



d 



z = cos 9 . 



Since T"'' = T'"', the contraction of {d/dz)i, with Eq. ^ leads to 



dz 



(35) 



(36) 



b) 



where the braces "( )" in the subscripts denote symmetrization of tensor. Equation 
governs the dynamical equilibrium in the 6'-direction. 
It is easy to check that 



V(a 



dz) 



(terms (x z) + Oi^z^) 



(37) 



b) 



near the equatorial plane {z = Q). The notation 0{x"') denotes terms of order or higher. 
Since T"'^ is regular and nonzero on the equatorial plane, from Eqs. (^^1) and (|n7|) we have 

d 



dz 



T 



ab 



oc z 



(38) 



for small \z\. 



For our model, T^' 



and = B'^ j'^ix, so we have 



d_ 

dz 



T 



ab 



Sir 



(39) 



where S is defined by Eq. ()3}. Then, by Eqs. ()H8j) and we have 

5^ = (terms independent of 6) + 0(cos^ 6) (40) 

for cos^ 6 ^ 1. Equation ()40|1 indicates that to the first order of cos^^, the dynamical 
equilibrium in the 6'-direction is guaranteed though the variable 6 does not appear in our 
solutions. Indeed, the solutions are exact on the equatorial plane 6 = 7r/2 . 



C. The radial momentum equation 

Using = g^tu*" + g^^u'^ and Ut = gttu^ + gt<t>u'^ , we can solve for and from Eqs. 
and (jHUj) . The results are 

t _ (/e - ^/l) + ^ (/e - ^^/l) 

^ i [^^/e + (1 - ^) /l] + ^^^^ (/e - ^^^/l) 
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Here and in some of the following equations we substitute the frame dragging frequency uj 
for the specific angular momentum a to simplify expressions. 
The corresponding covariant components are 

/e + ;i (/e - ^^k) (1 - ^ + ^^f^*) 



E 



(43) 
(44) 



The covariant radial velocity is Ur = QrrU^ = r'^u^ /A . 

Having expressed all components of and Ua in terms of u^, we are ready to write down 
the radial momentum equation for m*" by making use of UaU"' = — 1. The result is 



[ru 



-A 



(/e - n^k? 



^(uj~n^)k+{i-^ + ^^^^)k 



X f.2 



(45) 



The radial momentum equation (jl^j) is a nontrivial algebraic equation: it has multiple 
roots and singularities. Any physical solution must smoothly pass these singularities, which 
gives rise to strong constraints on physically allowable integral constants. As we will see in 
the next section, among the four integral constants appearing in Eq. (jUj), cq, /e, and 
/l, only three of them are independent for physical solutions. 



IV. ANALYSIS ON THE RADIAL MOMENTUM EQUATION 

In order to correctly solve the radial momentum equation, we must first understand the 
analytical features and the asymptotic behaviors of the solutions. In this section we analyse 
the singularities in the radial momentum equation (Sec. IIV A|) . and study the solutions on 
the horizon of the black hole (Sec. IIVB|) and at infinity (Sec. IIV C|) . 



A. Critical points 

Though [x' - (A/r^) {u - n^f] appears in the denominators on the right-hand side of 
Eq. (jl^ . it does not represent a singularity of the equation. Indeed, the factor disappears 
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from the denominators if we expand the second term on the right-hand side of Eq. ()45p . 
then combine with the first term. However, the factor 



1 + 



(46) 



which also appears in the denominators on the right-hand side of Eq. (jl^ . may represent a 
singularity of the equation. 



Introducing the relativistic Alfven four-velocity Cj^ = B"-/ \/ A^Txp^y + 5^, which satisfies 
cxaC^ < 1, we have 



1 + 



1 



1 - r2 



1 - 



UrW 



(47) 



Hence, the factor in Eq. ()4()|1 represents a singularity at the Alfven critical point defined by 

UrU^ — CArCA = , at r = Ta . (48) 



If we define a generation function 



F fr, u' 



1 + 



+ 



{UJ - /e + (1 - 2M + A^n^^ f^-^ ■ 



(49) 



then the radial momentum equation (j45p is given by F{r^u'^) = 0. The corresponding 
differential equation for can be derived from 

dF du'' OF 



From Eq. (jl^ . we have 



OF 



2 




















H 









du^ dr ~^ dr 



A 



(50) 



(51) 



Since 5^ = Anp^cl/ (1 - 4) and 3^3'' = Anp^CArC// (1 - 4), by Eqs. ^ and ^ 
we have 



j.2y^r 



A 



X 



[uJ - n^)^ - (/e - f2^/L) 



1 - r2 



(52) 
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Substituting Eqs. (jlTjl and into Eq. (jHU, we obtain 

dF 2 CKrC^ 



(1 - 4) 



1 



UrU 



UrU — 



r2 



(53) 



For a regular and stationary flow solution, cannot be zero at a finite radius since 
otherwise the mass density pm will diverge there. Therefore, Eq. ()53|) implies that there are 
two singularities in the differential radial momentum equation: one is at the Alfven critical 
point defined by Eq. pH|l . the other is at the fast critical point defined by 



UrU 



1 - r2 ' 



at r = rf 



(54) 



In a general magnetized wind theory, the Alfven critical point is not an X-type singularity 
and it does not inapqse any additional conditions on the solution for u'" except setting the 



integral constants 
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, 4^. For the model studied here, in fact the Alfven critical point is not 
a singularity at all, after the integral radial momentum equation (jl^ is taken into account. 
To see this, notice that 



dr 1 r'^u'^ 



+ 



— (u; - Vl^,) /e 



2M A 
h fi^ 



oc 



X 



— (CJ - ilq,) 



where in the last step Eq. (jl^j) has been substituted. Since the factor in Eq. appears 
in both dF/du'' and dF/dr, it disappears from the final form of the differential radial 
momentum equation (|5(J|) . Then, only the fast critical point, which is defined by Eq. (|54j) . 
is a true singularity in the radial momentum equation. 

Any physical solution must smoothly pass the singular fast critical point. So, by Eq. (jSHI), 
any physical solution must satisfy 



dF 



du^ 



, -:- 



r=r{ 



dF 
dr 







(55) 



r=r{ 



simultaneously. Equation (j55|) implies that the fast critical point corresponds to an ex- 
tremum of the generation function — indeed it is a saddle point of F as we will see in Sec. IVI 
Equation ()55j) together with F{r,u^) = allows us to solve for rf, u^ = u'^{r = rf), and 
/l for given cq, fl^^,, and /e. This implies that among the four integral constants Cq, fi^, 
/e, and /l only three are independent. 
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B. Solutions on the horizon of the black hole 



On the horizon of the black hole, where r = th, we have A = x = 0, = Qu, and 
A = AM'^r^, where = a/(2MrH) is the angular velocity of the horizon. Then, when 
7^ Qu, at r = th we can drop the first term on the right-hand side of Eq. and solve 
for My = u'^(r = th). The result is 

2M , „ „ _ /2M^ ^ 



= (/e - ^uh) + c'oi^u- , (56) 

where we have used the inequality /e — ^^h/l < 0, which comes from the requirement that 
on the horizon of the black hole energy must flow into the horizon locally. 

When Qxj, = Qn, it can be checked that as 5r = r — th — > we have oc 5r, uj — Qn oc 5r, 
1 - 2M/r + AuoVLH/r"^ oc 5r, and A/x^ = A/r"^ = AM'^ . Hence, when fivj, = ilu, on the 
horizon of the black hole we can drop the second term on the right-hand side of Eq. ()45p. 
Then we obtain u^^, which is just given by Eq. ()56|) with Qq, = fin- 
Substituting Eq. into Eqs. (jinj) and (HH), we obtain the specific energy and the 
specific angular momentum of fluid particles as they reach the horizon 

Eh = -/e + — cl (fin - , Lh = -/l + — (fin - f^*) . (57) 

From Eq. is always finite. However, UrU^ = (r^/A)(u'~)^ — *• oo as r — >■ r-^. Then, 

Eq. implies that B'^/Anp^ = c\/ (1 — c^) is finite at r = rn- So we have 

4/(1 - 4) 



UrU^ 



= . (58) 

Equation (j58p implies that fluid particles must supersonically fall into the black hole. 

The ratio of the electromagnetic energy density to the kinetic energy density, as measured 
by a LNRF observer, can be calculated by 



c 



87rp^r(r - 1) r(r - 



4 ^'^ 2 / o n2 



(59) 



where B'"^ and E'"- are given by Eqs. ()14|) and ()15|). F = —UaU^ is the Lorentz factor of fluid 
particles. As r — rn, we have F — > oo and 

cg(-.a)(nH-n.)^ (60) 

(/e-Oh/l)" 

where we have used (Fx)^^,.^^ = -E^h — -^^h^h = — /e + ^h/e • Hence, the ratio of the 
electromagnetic energy density to the kinetic energy density of the fluid is finite on the 
horizon of the black hole. 
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C. Solutions at infinity 



In Sec. IIII Al we have shown that when \1/ 7^ the solutions have a very bad dynamical 
behavior as r ^ 00: B'^/Airp^ oc {—u^')^^ 00, indicating that solutions with nonzero \l/ 
cannot extend to infinity. Therefore, here we set \1/ = to discuss the asymptotic behavior of 
the solutions at infinity. Then, from Eq. (jHUj) we have /e = —E. We assume that < E < 1, 
then we have /| < 1. 

When /I < 1 (i.e., 0<E <l), from Eq. (gH) we have F{r, m^) ^ (1 + cg/r^w^)^ (1 - /|) + 
fl/r'^ , as r ^ 00. Since the right-hand side is always positive, solutions to F = do not 
exist. This means that a fiow with /| < 1 cannot extend to infinity so must have an outer 
boundary of finite radius. 

When /I = 1 (i.e., E = 1), from Eq. (jl^ we have 



F{r,u'') ^ 1 



+ ^ ^ (61) 



as r 00. Then, there are two possible solutions to F = at infinity. One is ~ 
-(2M/r)i/2^C(r-3/2)^ ^j^j^l^ Eq. ^ leads to ru't' oc r'^ and oc {M/rf/'^ < 1. 

For the problem that here we are interested in, we expect that \ru'^ ^ 1 at large radii, 
so this asymptotic solution for is not what we are looking for. 
The other asymptotic solution for is 

+ , r^oo. (62) 

Then, from Eq. we have ru'^ ^ —{fh/r)/ (1 + c^/r'^u^) ^ [2M/r — {vJ'Yf^'^ as r — > 00, 
where in the last step Eq. and F = have been used. Then, the asymptotic ru'^ 
corresponding to Eq. (jU^ is 

/2M\ 

to"^ ^ ( j , r ^ 00 , (63) 

which leads to the correct asymptotic behavior |ru'^/u^'| oc (r/M)^/^ ^ 1. Note, Eq. ()63p 
corresponds to a super-Keplerian disk at large radii [a Keplerian disk has ru'^ = (M/r)^/^]. 
By Eq. ^,Eq. ^ implies that 

Pm ~ ^ ("^) ^ '^^^^^ ' ^ 00 • (64) 
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Since T ^ 1, T - 1 ^ (ru<^)V2, B'^ + E'^ ^ ^ {Cqu'*' /ru^'f as r oo, from Eqs. (jHl, 
dEH), and dSni) we have 

C ~ 1 , as r — i> oo . (65) 

Thus, at large radii the energy of the magnetic field is about equipartitioned with the kinetic 
energy of the flow. This explains the super-Keplerian nature of the solutions: at large radii 
the magnetic field is dynamically important, the "hoop stress" of the toroidal magnetic field 
provides an inward force confining the flow in addition to the gravitational force. 

The asymptotic solutions as r — * th and oo, for the case of \1/ = and E = 1, are 
summarized in Table IH 

V. SOLUTIONS WITH CONSTANT SPECIFIC ENERGY 

From the analysis in Sees. ITTl and IIIII we see that, for the model studied here, the integral 
constant has to be zero in order that the solutions can extend to infinity with finite r^B"^ 
and i?^/47rpm- For a flow with \E' = (i.e., = 0), the specific energy of fluid particles 
keeps constant, from r = oo down to r = rn. This type of solutions are of particular interest 
since they represent a new accretion mode in which accreting material loses its angular 
momentum but has its energy unchanged. 

From Sec. II VI at large radii the flow must be in the subsonic state since UrU^il — c^) / c\ ~ 
{M/rY <^ 1 . On the other hand, the flow must be in a supersonic state as it reaches the 
horizon of the black hole [Eq. (|58|)]. Therefore, a fast critical point [defined by Eq. (|54j)] 
must exist somewhere between r = rn and r = oo. Like in all wind problems, any physical 
solution must pass the fast critical point smoothly. 

In this section we solve the radial momentum equation to look for solutions corresponding 
to a flow accreting from infinity onto a black hole with constant specific energy. 

A. Solutions around a Schwarzschild black hole 

We study the simplest case first: solutions with constant specific energy around a 
Schwarzschild black hole (i.e., a = 0). The effects of the black hole spin on the solutions are 
studied in the next subsection. 



18 



When a = \1/ = 0, the generation function [Eq. ()49p ] becomes 



F fr, 



1 + 



2M 



Je 



1 _ 2M 



f2 



(66) 



The corresponding derivatives are 



dF 

dF 

dr 



w 

2 

r 



1 + 



2d 



2M 

r 

2M" 



X _ 2M 



2M 



(67) 



r 



1 - 



2M 



5M 



j.2r^r 



3M 



2/^ 



(68) 



As shown in Sec. IIV A| the only singular point in the radial momentum equation is the 
fast critical point that is defined by 



\ _ 2M 



y,2yr 



r 







from which we can solve for the radial velocity at r = rf 



r „2 



1 



2M 



. 2M , 
-1 + + /] 



r = Tf 



1/3 



(69) 



(70) 



Since m[ must be negative and rf > rn = 2M, the following condition must be satisfied 

2M 



- 1 



rf 



+ /e>0 



(71) 



Note that m[ does not depend on /l- 

From the integral radial momentum equation F{r,u^) = 0, we can solve for /l 



/l = -r 



1 + 



2M 



I — 2M 



- 1 



1/2 



(72) 



Substituting Eq. ^ into Eq. 



we can eliminate /l from Eq. 



and obtain 



dF 
dr 



1 + 



iu' 



j,2y^r 



/e 



1 - 



2M 

r 

3M 
r 



1 - 



j.2y^r 



y,2^r 



AM 

r 

2M 

r 



(73) 



For /l to be real, the term in the second pair of brackets on the right-hand side of Eq. ()72p 
must be nonnegative 



/I 



X _ 2M 



1 > . 



(74) 
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Making use of Eqs. (fTOj) and (f7T|) . at r = rf Eq. (fflj) becomes 

Equation (f75j) restricts the region in the (rf , cq) space where physical solutions exist. 

At the fast critical point Eq. must be satisfied. From Eq. (jU7j) . at r = rf dF/du'' = 
leads to Eq. ()7()|1 . From Eq. (ffH|) . at r = rf dF/dr = is equivalent to 



4 f AM\ , {u\Y - fi f 3M 



+ 



TfMf V J (] _ V '^f 



eg 2M 



rfu^j \ rf 



. (76) 



From Eqs. (fTOj) and (fTBj) we can numerically solve for rf = rf (/e, Cq) and Mf = m[ (/e, Cq) 
for any given /e and Cq. The location of the fast critical point in the (r, u"^) phase space is 
then determined. Substituting the solutions into Eq. (ff^ . we can determine /l = /l (/e, cq), 
which has a negative sign. 

From Eq. (jHUj) . we have /e = — -E when \l/ = 0. When the flow has an outer boundary 
with a finite radius tq, we assume that the constant specific energy of fluid particles, E, is 
equal to that of a test particle moving on a Keplerian circular orbit at r = rg. Then, we 
have /e = — -EKepierian('"o)- When ro = oo, we assume that /e = —E = —1. 

In Kg. rawe show the .olutio,. fa _ndi„g to ^ , where ^ 6 . the 
radius at the marginally stable circular orbit rg = M is the gravitational radius of the 
black hole. This represents the solutions in the plunging region from the inner boundary 
of a thin Keplerian disk to the horizon of the black hole. The solutions for rf, represented 
by the solid curves in the figure, consist of three branches: A, B, and C, intersecting at 
the dark point. The dashed curve represents the boundary for physical solutions to exist: 
above the dashed curve Eq. (ffH|) is violated so physical solutions do not exist. (The dashed 
curve is indeed coincident with the solid curves B and C. To make it visible, in the figure 
we have shifted the dashed curve upward a little bit.) From Fig. ^we see that, below the 
dark point (where r = 3.4477 rg, cq = 4.4030 rg), rf is uniquely determined by the value of 
Cq. As Cp/rg — > , rf approaches the radius at the marginally stable circular orbit r^s (the 
vertical dotted line). This is just what we expect for a weakly magnetized thin Keplerian 
disk: the (magneto) sonic point should not be far from the marginally stable circular orbit 
(13, 44,0] and references therein). Above the dark point, multiple solutions for rf exist. 
In this paper we focus on the situation of weak magnetic fields with Cg/rg ^ 1, so we do 
not discuss the solutions above the dark point in detail. 
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We are more interested in solutions with tq ^ rg, corresponding to a magnetized flow 
accreting onto a black hole from large radii. In Fig. |21 we show the solutions for rf corre- 
sponding to To = lO^rg (dark lines) and Tq = oo (gray lines). The solutions for Tq ^ rg have 
the same topology as that for tq = r^^- However, when Co/vg is fixed, rf decreases as 
increases. For ro ^ rg, rf becomes insensitive to the value of tq: on branches A and B the 
solutions for rg = lO^rg are almost indistinguishable from those for ro = oo. This is caused 
by the fact that /e ~ —1 for ro ^ rg. From Fig. |21 we see that, for the case of ro ^ rg, 
rf approaches the radius at the marginally bound circular orbit r^b = 4rg |^ (the vertical 
dotted line) in the limit cl/r"^ . 

In both Figs. ^ and |21 the dashed curve is coincident with the solution branches B and 
C. So, the solution branches B and C correspond to /l = [see Eqs. (f7^ and (f7i|) ]. By 
Eq. (jHZj), such solutions correspond to Lh = since = Qu = 0. Thus, for a sufficiently 
strong magnetic field (corresponding to the branches B and C), fiuid particles lose all of 
their angular momentum as they reach the horizon of the black hole. 

Once Tf, Mf, and /l are determined as functions of cq and /e, we can solve the radial 
momentum equation for the radial velocity u^' as a function of radius r. However, generally 
it is not a simple thing to solve the algebraic equation F{r, u"') = 0, which has singularities 
and multiple roots. In practice we would rather solve the first-order differential equation 

1-^(1- 2M\" 
r'^u^ \ r J 

^1^^ " iurf (1 - M£)- _ _iL (1 _ _ /I) ' 

(77) 

which is obtained by substituting Eqs. (jU7|) and (fTSj) into Eq. fl50|) . Note that the only 
singularity in Eq. (|77j) is the fast critical point. 

We choose on each side of the fast critical point a pair of (rj, «[), which can be obtained 
by solving the algebraic equation F{r, u^) = where initial guessed values can be obtained 
by visually checking the contour plot of F versus r and u^. Then, starting from the point 
(rj, u[), we integrate Eq. (f77j) inward and outward, toward the horizon of the black hole and 
the fast critical point if (rj,M[) is on the left-hand side of the fast critical point, or toward 
the fast critical point and large radii if (rj,-u[) is on the right-hand side of the fast critical 
point. [The starting point (r,, M[)cannot be exactly at the fast critical point since the latter 
is a fixed point of the differential equation.] 



dlnu^ 1 - ;i (1 - ^) + [{un' - M] (1 - '-^r (1 - ^) 
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As an example, we have solved the equation with /e = —1 (i.e., ro = oo) and cq = 0.1 rg, 
corresponding to rf = 3.98852 rg, m[ = —0.053964, and /l = —3. 96522 rg. The contours of 
the corresponding generation function are plotted in Fig. El I59I], from F = —0.2 to F = 0.2 
with step 6F = 0.05. The contours of F = 0, which correspond to the physical solutions, 
are the two diagonal lines that smoothly pass the fast critical point at the dark point. The 
solution that we are interested in here should have < —u^ ^ 1 at large radii and ■u'' ~ — 1 
near r = th, which corresponds to the diagonal line that goes from the right-bottom corner 
to the left-top corner in Fig. El The differential radial momentum equation ()77|) is then 
solved along this diagonal line, with the approach described above. 

The solution for the radial velocity function of r is plotted in Fig. HI (upper 

panel). For comparison, we show the asymptotic solution as r — cxd with the dashed line 
[see Eq. flH^ ]. We have calculated C, = [urU^' (1 — c\) /c\]^^'^, and shown the results in the 
lower panel of Fig. HI The parameter ^ measures the ratio of the radial velocity to the Alfven 
velocity. At the fast critical point we have ^ = 1 [Eq. (j5l|) ]. The figure clearly shows that 
the flow starts from a subsonic state < 1) at large r, passes the fast critical point at r = rf 
(^ = 1), then enters a supersonic state (^ > 1) at r < rf. The ratio ^ blows up at r = rn, 
confirming our analytic analysis in Sec. IIVBI [see Eq. (j^Hj) ]. 

The specific angular momentum of the fiuid particles, L, can be calculated by substituting 
a = \1/ = 0, and the solution u^{r) into Eq. ()44|) . The results are shown in Fig. For 
comparison, we show the specific angular momentum of a thin Keplerian disk with the 
dashed line. The Keplerian disk is joined to a free-fall fiow in the plunging region, so its 
specific angular momentum keeps constant inside the marginally stable circular orbit. The 
fast critical point and the marginally stable circular orbit are indicated by two arrows (left 
and right, respectively). From Fig. El we see that, a disk fiow with constant specific energy 
can lose angular momentum at a rate similar to that of a thin Keplerian disk. 

Figure El shows that the specific angular momentum of the fiow remains almost constant 
after the fiow passes the fast critical point. This suggests that we can treat the fast critical 
point as an inner boundary of the accretion fiow. To justify this, we have calculated the mass 
density of the fiow and the ratio of the radial velocity to the rotational velocity, as functions 
of radius r. The mass density of the fiow is simply given by Eq. ()27|). which approaches a 
constant as r ^ 00 [Eq. ()64|) ]. The ratio of the radial velocity to the rotational velocity, as 
measured by a local static observer, is given by Vr/v^ = ru^/xL . The numerical results for 
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the mass density and the ratio of the radial velocity to the rotational velocity, corresponding 
to the solution in Fig. HI are shown in Fig. IHl 

We see that, at large radii |fr-/f0| <^ 1, and pm ~^ Poo = F"^/ {^'^CD [Eq. As 
the flow gets close to and passes the fast critical point (the dark points in the figure), 
Pm/ Poo drops quickly and \vr/v^\ increases quickly. At the fast critical point we have Pm = 
0.01167poo and Vr = — 0.07629 f<^. These numbers support the point that the fast critical 
point behaves like an inner boundary of the disk accretion flow: it marks the transition of the 
flow from a state with \vr/v^\ <^ 1 at large radii to a state with \vr/V(f,\ > 1 near the horizon 
of the black hole. Beyond the inner boundary, the magnetic field plays an important role 
in the dynamics of the flow and transports angular momentum outward. Inside the inner 
boundary, the dynamical effects of the magnetic field become unimportant and the fluid 
particles free-fall to the horizon of the black hole (Fig. E}. 

Figure also shows that, as r — th, Pm is finite but \vr/vtf,\ blows up. The former is 
because of the fact that u"^ is finite at r = rn (see Sec. IIVB|1 . the latter is caused by the 
well-known fact that and ^ — 1 as particles approach the horizon. 

The ratio of the magnetic energy density to the kinetic energy density of the flow, as 
measured by a local static observer [see Eq. (jSHI); note that when a = \E' = the electric 
field is zero, see Eq. ()15p]. is calculated and shown in Fig. [7| In the region of r ^ rf, the 
ratio C is ~ 1, i.e., the energy of the magnetic field is about equipartitioned with the kinetic 
energy of the flow [Eq. ()65|l ]. This is caused by the fact that for r ^ rf the velocity of the 
flow is predominantly rotational: |fr/f<^| -C 1 (Fig. IH)), so the magnetic field is sufficiently 
amplified by the shear rotation until the energy density of the magnetic field is comparable 
to the kinetic energy density of the flow. In the region of r < rf, we have C 1, caused by 
the fact that for r < rf the radial component of the flow velocity becomes comparable to, or 
even greater than the azimuthal component (Fig. El), then the amplification of the magnetic 
field by the shear motion becomes unimportant [35]. At the fast critical point (the dark 
point in the figure), we have ( = 0.004954 . As r — rn, we have ( . 

To see how sensitive the solutions are to the variation of cq, and how the solutions 
approach the asymptotic solutions at infinity found in Sec. IIV CI we have solved the radial 
momentum equation with /e = —1 corresponding to two different values of Cq: Cq = Irg 
and Cq = 10~^rg. The solutions for the angular velocity Q = u'^/u^ are shown in Fig. |H1 
The solid line corresponds to cq = 1 rg, the dashed line to cq = 10"^rg. The dotted line in 
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the figure is the angular velocity of a thin Keplerian disk joined to a free-fall flow in the 
plunging region. The dashed-dotted line shows the asymptotic solution for Q derived from 
Eq. (jnSl)- We see that, the solutions are insensitive to the value of Cq for Cq < r^, and they 
are well represented by the asymptotic solution in the region with log(r/rg) > 1.5 . 



B. Solutions around a Kerr black hole 



The spin of the black hole has important effects on the solutions of an accretion flow near 
the horizon of the black hole. For example, for the case of a thin Keplerian disk, the spin of 
the black hole determines the location of the inner boundary of the disk, how much energy 
is dissipated in the disk, and how much angular momentum is carried into the black hole 
by the flow [3, 0, 0]. In this subsection we show how the spin of the black hole affects the 
accretion process by presenting some solutions for a magnetized flow with constant speciflc 
energy accreting onto a Kerr black hole. We assume that at large radii the speciflc angular 
momentum of the disk is always positive, but the speciflc spinning angular momentum of 
the black hole, a, can be either positive or negative. 

For the case of a Kerr black hole with \1/ = 0, the generation function becomes 
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The derivative of F with respect to is 
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(79) 



The derivative of F with respect to r is too lengthy so we do not write it out here. 



The radial velocity at the fast critical point can be solved from dF/ du^ 
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From the radial momentum equation F{r,u^) = 0, we can solve for /l 
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Substituting Eq. ()81|) into the expression for dF/dr, we can eliminate /l from OF /dr. Then, 
similar to the case of a Schwarzschild black hole, we can obtain rf and u[ by solving the 
joined equations (jSUI) and dF/dr = 0. Then by Eq. (|HT|) we can determine /l = /l(co, /e)- 
For physical solutions /l must be real, so the following condition must be satisfied 

?^^),™.f-A(l-a^-/l),o. (S2, 

Making use of Eq. (jHUj) , at r = rf Eq. (jH^ becomes 

+ /^ • (83) 
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As in the case of a Schwarzschild black hole, Eq. (jHH|) provides a constraint on the region in 
the (rf , Co) space where physical solutions exist, whence /e is specified. 

In Fig. El we show the solutions for rf with ro = rins(a.) and /e = — -EKopiorian('"05 ct); where 
EKepieTia.n{ro,a) is the specific energy of a test particle moving on a Keplerian circular orbit 
around a Kerr black hole at r = rg 0, J, |^ . Each panel corresponds to a different spinning 
state of the black hole. The dashed curve shows the boundary for physical solutions: above 
the dashed curve Eq. (jEHI) is violated so physical solutions do not exist. For the cases of 
a/M = 0.95 and 0.998, we have rms < 2M and Eq. (jHH|l is always satisfied since the right- 
hand side is always positive [see Eq. (jHOI) where m[ < 0]. So, the dashed curve does not 
appear in the panels corresponding to a/M = 0.95 and 0.998. In the case of a/M = —0.9, 
there are two branches of solutions for rf, labeled by A and B. This is a general feature for 
the case of a < 0. Like in the case of a Schwarzschild black hole, in the limit of Cg/rg 0, 
rf approaches rms, except for the case of a/M very close to 1, e.g., a/M = 0.998. In the 
case of a/M = 0.998, solutions for rf do not exist if c^/r"^ is sufficiently small. This means 
that for sufficiently small Cg/rg, the flow is already supersonic as it leaves the marginally 
stable circular orbit. This is caused by the fact that we have neglected the gas pressure in 
our treatment. In the Keplerian disk region (r > rms) the gas pressure cannot always be 
neglected — especially near the inner boundary (r ^ r^^). When the gas pressure is included, 
it will make the radial motion subsonic in the disk region and give rise to a sonic point close 
to the marginally stable circular orbit 

(0 

In Fig.^Jwe show the solutions for rf with rg = oo and /e = —1. Similar to the case of a 
Schwarzschild black hole, when c^/r^ -C 1 the fast critical point approaches the marginally 
bound circular orbit: rf ^ r^ih, where rmb = 2M [l — a/2M + (1 — a/M)-*^/^] ^S,]. 
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By comparing Figs. M and HUl with Figs. [T] and we can see how the topology of the 
solutions for rf evolves with the spin of the black hole. For a < 0, the solutions have two 
branches: A and B. As a/M increases from —1 to 0, the branch B moves toward left. (One 
can check this by solving rf for a negative a > —0.9.) As a approaches 0~, the branch B of 
negative a disappears, the branch A of negative a becomes the branches A and C of a = 0. 
For a > 0, the solutions have only one branch. As a/M decreases from 1 to 0^, the branch 
of positive a becomes the branches A and B of a = 0. 

As examples, we have solved the radial momentum equation for a Kerr black hole with 
a/M = 0.95 and a/M = —0.9, respectively. We assume that cq = 0.1 rg and /e = —1 (i.e., 
tq = oo). In Fig. ^2 we show the specific angular momentum of fluid particles (upper panel), 
and the ratio of the electromagnetic energy density to the kinetic energy density as measured 
by a LNRF observer [lower panel, defined by Eq. ()59|) ] for the case of a/M = 0.95. In the 
upper panel, the dashed curve shows the corresponding solution for a thin Keplerian disk 
joined to a free-fall flow inside the marginally stable circular orbit. The two arrows show 
the positions of the fast critical point (left) and the marginally stable circular orbit (right). 
In the lower panel, in addition to the total ( (solid line), we also show the contribution from 
the magnetic field (dashed line) and the electric field (dotted line) separately. The dark 
point on the solid line shows the location of the fast critical point. The two circles show 
the location of the black hole horizon. Solutions for the specific angular momentum and 
the ratio of the electromagnetic energy density to the kinetic energy density for the case of 
a/M = —0.9 are shown in Fig. IT^ 

From Figs. ^2 and El we see that, similar to the case of a Schwarzschild black hole, the 
specific angular momentum of the fluid particles is effectively removed by the magnetic field, 
though the specific energy keeps constant during the motion. However, the specific angular 
momentum of the fluid particles as they reach the horizon of the black hole depends on the 
spin of the black hole, which can be seen by comparing Figs. ^2 and El with Fig. El 

At large radii the energy of the electromagnetic field is about equipartitioned with the 
kinetic energy of the particles, while near and inside the fast critical point the ratio of the 
electromagnetic energy density to the kinetic energy density drops to values ^ 1. However, 
unlike the case for a Schwarzschild black hole, the ratio is not zero on the horizon of the 
black hole when a 7^ 0. This is related to the fact that when a 7^ the magnetic field is 
not aligned with the velocity field near the horizon (caused by the frame dragging effect. 
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see Sec. |n}, resulting a nonzero electric field as shown by the dotted lines in the two figures 
(lower panels). On the horizon the ratio of the electric energy density to the magnetic energy 
density approaches 1, so the dashed line and the dotted line end on the same small circle. 
The strength of the electric field drops quickly as the radius increases, and at large radii 
(where E'^E'"" oc r^*^, B'^B"^ oc r^^) the effect of the electric field is negligible. 

The fact that the ratio of the electromagnetic energy density to the kinetic energy density 
is small near and inside the fast critical point, for both the Schwarzschild case and the Kerr 
case, explains why the fast critical point approaches the marginally bound circular orbit in 
the limit of tq ^ and Cg/rg <^ 1, because then E 1 and the dynamical effects of the 
electromagnetic field are unimportant near and inside the fast critical point. The results 
support the point that the fast critical point behaves like an inner boundary of the fiow, 
beyond which the dynamical effects of the magnetic field are important, inside which the 
dynamical effects of the magnetic field are negligible and fiuid particles free-fall to the black 
hole. 

In Sec. lIVBI we have shown that the specific angular momentum of the fiuid particles on 
the horizon of the black hole is given by Eq. ()57|1 . For solutions that smoothly pass the fast 
critical point, the parameter /l is determined as a function of /e and Cq when \E' = 0. Hence, 
if /e is specified, for a given spin of the black hole Lh depends only on cq. In Fig. we 
show Lh as a function of Cq, where different curves correspond to different spinning states of 
the black hole. The outer boundary of the flow is assumed to be at inflnity so that /e = —1. 
Figure Uni conflrms the results in Figs. El CH and[T21 Furthermore, we see that Lh depends 
on a in opposite ways for small and large cq: for small cq, Lh decreases with increasing a; for 
large Cq, Lh increases with increasing a. And, for sufficiently large Cq, Lh becomes zero when 
a = (see Sec. IV A|) . For large cq and negative a, Lh is negative. The Lh corresponding to 
positive a is always positive. 

VI. DISCUSSION: THE EFFECT OF FINITE RESISTIVITY AND DIFFERENT 
EVOLUTION OF LARGE AND SMALL-SCALE MAGNETIC FIELDS 

The solutions that we have found have the following unique feature: the magnetic fleld 
is aligned with the velocity fleld [Eq. and, as a result, at large radii the energy of the 

magnetic fleld is equipartitioned with the kinetic energy of the flow [Eq. ()65|)]. The latter is 
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obvious from the following consideration: At large radii, where the gravity is approximately 
Newtonian and the radial flux of angular momentum is not important, we have AtcPt^u^u'^ ^ 
B^'B'^. Then, B''/B^ = ju"^ implies that B^B"" ^ ^ixp^UrU^ and B^B-^ ^ ^■np^u^u'^, 
i.e., ^VStt ^ Pmv'^/'2,, where B^ = B,.B'' + B^B'^ and v"^ = Ur-u'' + u^u'^. An almost 
inverse statement is also true: If at large radii B^B"^ At^ p^u^u^ , then we must have 
B^/B'^ ^ u'^/u'^ at large radii. (See App. iBlfor the corresponding solutions for a Newtonian 
accretion disk.) 

The above results rely on the fact that we have neglected all energy dissipation processes 
by assuming a perfectly conducting fluid with zero resistivity and viscosity. In the Newtonian 
limit (which is valid at large distances from the central black hole), the induction equation 
including a finite electric resistivity rj is given by (see, e.g., jlS^ ) 

^=Vx fvxB-^VxB) . (84) 



dt V 47r 

We assume that r] = constant. Then, substitute Eq. ()B5|) into Eq. (jHD), we obtain the 
stationary and axisymmetric solution 

r {vrB^ - v^Br) = ^4- (^^*) + C , (85) 
47r ar 

where C is an integral constant. 

To estimate the effect of finite resistivity, let us consider the case of C = 0. Then from 
Eq. we have 

Br ^ Vr ^ n din irBs) , s 

7^ = Cm-, Cm = 1--^ ir^- 86 

B^ V(f, AirrVr amr 

The parameter (u determines the orientation of the magnetic field relative to the velocity 
field in the r — plane. Define 

_ 87i^r\v^\ _ dlnjrB^) 

-Km = , q= — Tj , [oi) 

rj alnr 

where Rm is the magnetic Reynolds number associated with the azimuthal motion, q mea- 
sures the radial gradient of the magnetic field. Then, (u can then be rewritten as 



Vr 



(88) 



A 2nq 

Therefore, the effect of finite resistivity is determined by the parameter A. When |A| ^ 1, 
we have Cm ~ 1, the effect of resistivity is unimportant and the solutions correspond to 
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a perfectly conducting fluid as in the case of our analytical model. Then, the magnetic 
field is aligned with the velocity field and the magnetic energy density is equipartitioned 
with the rotational energy density of the disk at large radii. When |A| ^ 1, we have 
ICmI ~ l/l-^l ^ 1) the effect of resistivity becomes important and the solutions correspond 
to a nonperfectly conducting fluid. The magnetic field is not aligned with the velocity field, 
resulting significant energy transportation and dissipation. Then, in the stationary state the 
magnetic energy density should be much smaller than the rotational energy density of the 
disk (as shown bellow). The effect of finite resistivity is particular important for small-scale 
and chaotic magnetic fields that reverse directions frequently, since at places where = 
but dB^jdr 7^ we have |g| = 00 [Eq. (jHTj)] and A = 0. 

It is the assumption of a perfectly conducting fluid with zero resistivity and viscosity that 
makes our analytical results different from the MHD numerical simulation results, where 
the magnetic field is not always aligned with the velocity field and the energy density of the 
magnetic f ield is usually much smaller than the rotational energy density of disk particles 
22, |2^ Q, In numerical simulations, finite resistivity and energy 



(see, e.g., |2JJ 



dissipation arising from finite spatial gridding in numerical codes always exist and are o 
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ten 

m. 



large, even though the physical resistivity and viscosity can be assumed to be zero 
In fact, in current three-dimensional MHD simulations, the Reynolds number arising from 
finite numerical gridding is often smaller than the magnetic Reynolds number associated 
with most astrophysical plasmas 0, li^ ■ 

How about the effect of finite resistivity in real astronomy? For astrophysical plasmas, the 
magnetic Reynolds number Rm is usually extremely large (> 10^) due to small microscopic 
Ohmic resistivity and large spatial scales and velocities (see, e.g., {4^). However, dissipation 
processes like convection, turbulence, etc, can induce a large effective resistivity jsy, |51 , 



53|, thus produce a smaller magnetic Reynolds number. The parameter q can also be large 



if the magnetic field has small scales and is chaotic. 

For an ordered large-scale magnetic field we usually have g ~ 1 in regions where 7^ 
0. While for small-scale and chaotic magnetic fields, the parameter q can be very large, 
especially at places where B^ = 0. In this latter case, the accretion disk cannot always 
be described with a simple stationary and axisymmetric model without using a statistical 
approach. We can define a coherent length scale Ic so that |g| r/l^. Then, A in Eq. 
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can be recast as 



, Rc 

A ~ 

27r 



Vr 



K^'-^, (89) 

V 



where Vc is the velocity associated with the local coherent structure, Rc is the corresponding 
magnetic Reynolds number, and r] is the resistivity in the disk. 

Thus, accretion disks driven by large-scale magnetic fields behave very differently from 
those driven by small-scale magnetic fields. For an accretion disk driven by a large-scale 
magnetic field, like the model studied in this paper, A is usually ^ 1 (since g ~ 1) and 
determined by the global motion of the disk through Rm and Vr/vtf, [Eq. (jHH|) ]. For an 
accretion disk driven by small-scale magnetic fields, A is usually <^ 1 (since |g| ^ 1) and 
determined by Rc and Vr/vc [Eq. (jH^ ]. In this latter case, the global motion affects A (so 
affects the equilibrium state) only through the radial velocity Vr — the global rotation of the 
disk has no effect on A (so does not affect the equilibrium state of the magnetic field). 

Making use of Eq. (jHBj) and Anp^VrV^ ~ BrB^, we can estimate the energy density of the 
magnetic field by 

_ Cm +CM{Vr/Vc^? /Q^v 

where Su = (l/Svr) [Bf. + S^) is the magnetic energy density, Sk = (Pm/2) [v^ + f^) is the 
kinetic energy density of the disk. When A ^ 1 (corresponding to large-scale magnetic 
fields), we have Su ~ ^k, the magnetic energy is approximately equipartitioned with the 
rotational energy of the disk (recall that jfr/f^l ^ 1). When A ^ 1 (corresponding to 
small-scale magnetic fields), we have Sm/Sk ~ (-Rc/Stt) \vr/v^\ ~ A \vc/v^\ <^ 1, the mag- 
netic energy density is much smaller than the rotational energy density of the disk. As an 
example for an accretion disk with small-scale magnetic fields, consider a thin Keplerian 
disk with ~ v"^/ \vrh\ where Vc corresponds to the fiuctuation in the fiuid velocity arising 
from turbulence [1^, and -Rc ~ Stt (see, e.g., js^), then we have Su ~ Prn'^c/^) i-e., the 
magnetic energy density is approximately equipartitioned with the kinetic energy density 
in turbulence. Thus, for a disk with small-scale magnetic fields the equilibrium state is 
determined by the local turbulent motion rather than the global rotation [si. ^|. 

Based on the above arguments, we believe that the results of Gammie 2^ more 
relevant to a large-scale magnetic field connecting the disk to the black hole, rather than 
small-scale and chaotic magnetic fields in the plunging region. 

30 



VII. SUMMARY AND CONCLUSIONS 



With an analytical model we have investigated the dynamics of a cold accretion disk 
containing a large-scale magnetic field around a Kerr black hole, trying to understand the 
process of angular momentum and energy transportation. A one- dimensional radial momen- 
tum equation is derived near the equatorial plane of the black hole, which has one intrinsic 
singularity at the fast critical point. Any physical solution corresponding to an accretion 
fiow starting from infinity subsonically must smoothly pass the fast critical point in order 
to reach the horizon of the black hole. Because of this requirement, among the four integral 
constants (cq, fi^, /l, and /e) only three are independent: one of them has to be determined 
function of the other three. 

In order for r'^B^' and i?^/47rpm to be finite at infinity, the constant fi^r has to be zero 
(Sees. Inland IIlIA|) . This indicates that in the stationary state the magnetic field and the 
velocity field of the fiuid must be parallel to each other. When this condition is satisfied 
{Q^ = 0), the specific energy of the fiuid particles remains constant during the motion, but 
the specific angular momentum changes. 

The radial momentum equation has been solved with the assumptions that fi^, = 0, 
Cg/rg <^ 1, and /e = —1 (corresponding to an accretion fiow with the outer boundary at 
infinity). The results can be summarized as follows: 

1. The accretion fiow starts from a subsonic state at large radii, passes the fast critical 
point near the marginally bound circular orbit, then enters the black hole superson- 
ically (see, e.g.. Figs. EJ. The solutions are super-Keplerian (Figs. El HH and |T21 
upper panels), whose asymptotic behaviors on the black hole horizon and at infinity 
are summarized in Table HI 

2. The specific energy of the fiuid particles remains constant, but the specific angular mo- 
mentum is effectively removed by the magnetic field and transported outward (Figs. 
^2 and ^1 upper panels). The rate of transporting angular momentum is similar to 
that for a thin Keplerian disk, though the specific angular momentum that the fiow 
finally carries into the black hole is somewhat larger. 

3. At large radii the electromagnetic energy density is about equipartitioned with the 
kinetic energy density of the fiuid particles, but near and inside the fast critical point 
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the ratio of the electromagnetic energy density to the kinetic energy density is <^ 1 
f Figs. 171 im and IT^ lower panels). This is due to the fact that at large radii the velocity 
of the fluid is predominantly toroidal so the magnetic field is sufficiently amplified by 
the shear rotation, while near and inside the fast critical point the radial velocity 
of the fluid becomes important then the magnetic field has no time to get sufficient 
amplification (see Fig. IH)). 

4. When the black hole is Schwarzschild, the electric field measured by a local static 
observer is zero and the ratio of the magnetic energy density to the kinetic energy 
density approaches zero on the horizon of the black hole (Fig. [7j). 

5. When the black hole is Kerr, due to the frame dragging effect the electric field measured 
by a LNRF observer is nonzero. The electric field decays rapidly with increasing radius, 
but near the horizon the electric field makes a significant contribution to the total 
electromagnetic energy density. On the horizon, the ratio of the total electromagnetic 
energy density to the kinetic energy density, which is contributed equally by electric 
field and magnetic field, is tiny but nonzero (see Figs. [TT] and [T21 lower panels). 

6. The specific angular momentum of the fluid particles as they reach the horizon, Lh, 
is a strong function of a/M and Co/rg, though in the limit Cg/rg ^ 1 Lh becomes 
insensitive to the change in Cq/t^ (Fig. IT^ . For a positive a, Lh is always positive. 
For a zero or negative a, Lh becomes zero or negative for sufficiently large Cq/t^. 

With the existence of solutions with constant specific energy for a magnetized accretion 
flow, the question posed in the Introduction {Do magnetic fields transport or dissipate energy 
as efficiently as they transport angular momentum?) is at least partly answered: Magnetic 
fields do not have to transport or dissipate a lot of energy as they efficiently transport angular 
momentum. This has important implications to the theory of accretion disk: an accretion 
disk driven by magnetic fields may have a very low radiation efficiency. People generally 
believe that an accretion disk usually has a radiation efficiency much higher than a spherical 
accretion flow. The belief comes from the following argument: a disk has a lot of angular 
momentum which must be removed in order for the disk material to be able to accrete onto 
the central compact object. In the meantime, a lot of energy is expected to be removed 
also. The results in this paper show that the specific energy of fluid particles does not 
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have to change when its angular momentum is being removed. In fact the solutions with 
constant specific energy are the unique solutions for the model studied in this paper, all other 
solutions have bad asymptotic behaviors at infinity (namely, r^B^ and B'^/Airp^ diverge), if 
the specific angular momentum of the fiuid is unbounded at infinity. 

We have ignored energy dissipation by assuming that the disk plasma is perfectly con- 
ducting. Then, inevitably, at large radii the magnetic energy density is equipartitioned with 
the kinetic energy density, indicating that the accretion fiow must be geometrically thick. 
The solutions presented in the text apply only to the region near the central plane of the 
disk. However, in App. El we show that the solutions can be easily generalized to more 
general and more complicated models, including the disk region well above or well below 
the equatorial plane. So, solutions with constant specific energy exist for more general and 
more complicated models. 

The effects of finite resistivity are discussed in Sec. I VII When the resistivity is nonzero, an 
accretion disk driven by a large-scale and ordered magnetic field behaves very differently from 
that driven by small-scale and chaotic magnetic fields. For a large-scale and ordered magnetic 
field, in a stationary and axisymmetric state the magnetic energy density is approximately 
equipartitioned with the rotational energy density at large radii if the resistivity is small. 
For small-scale and chaotic magnetic fields that reverse directions frequently, even if the 
resistivity is small, in the equilibrium state the magnetic energy density is likely to be 
equipartitioned with the kinetic energy density associated with local random motions (e.g., 
turbulence). 

A finite resistivity in the disk will give rise to a nonzero radiation efficiency. However, a 
nonzero radiation efficiency may also arise in the following way: a disk corresponding to our 
solutions has a finite outer radius tq and beyond which it is joined to a thin Keplerian disk. 
Then, the arguments in Sec. IIII Al indicate that the region inside tq might have a radiation 
efficiency < ^qLq ~ Tg/ro , assuming that ro ^ rg . The thin Keplerian disk outside rg 
would give rise to a radiation efficiency ^out ~ ^g/i^o ■ Then, the total radiation efficiency of 
the disk region from r = rn to r = Too would be etotai = £m + £out ~ f^g/f^o ■ 

Finally, we remark that the issue of stability/instability has been ignored in our treatment. 
MHD instabilities are often important in the study of magnetized accretion disks (for reviews 



see 



IS 



53l|). they may also be important for our solutions. However, a detail treatment 



of the stability/instability of our solutions is beyond the scope of the current paper, so we 
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would leave it for future work. 
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APPENDIX A: GENERALIZATION OF THE SOLUTIONS TO MORE GEN- 
ERAL MODELS 

The results presented in the text can be generalized to more general models where the 
stream lines of the fluid do not lie on surfaces of constant 6. To do so, let us assume that 
the disk can be foliated with a set of surfaces so that each stream line of the flow lies on one 
of such surfaces. We call the surfaces so defined the stream foliation surfaces. By symmetry, 
the equatorial plane is one of the stream foliation surfaces. The upper and lower faces of 
the disk are two stream foliation surfaces. 

For a stationary and axisymmetric disk, the stream foliation surfaces are specified by some 
functions of r and 9. We can define a congruence of lines that perpendicularly intersect each 
stream foliation surface. Each point on the equatorial plane has a radial coordinate r, the 
congruence of lines carry this radial coordinate to each of the other stream foliation surfaces. 
We label each stream foliation surface with a variable G, and choose B = on the equatorial 
plane, 9 = ±1 on the upper and lower faces of the disk. Then, we are able to define two 
new coordinates R and G, which are related to r and 9 by 

R = R{r,9), e = e{r,9), (Al) 

so that each stream foliation surface corresponds to O = constant and the coordinate lines 
of R and B are orthogonal to each other. On the equatorial plane G = we have R = r. 
The relations in Eq. ()A1|) should give rise to a one-to-one correspondence between (r, 9) and 
{R, 6) at least locally. 
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A general stationary and axisymmetric metric can be written as 

ds^ = Qtt + 2gt^ dtd(f) + g^^ d(f + g„ dr'^ + geg d9^ , (A2) 

where gf^i, = g^^ are functions of r and 9. Now let us make coordinate transformations from 
(t, r, 6*, 0) to (t, i?, 6, 0). Then the metric becomes 

ds^ = gtt dt^ + 2gt^ dtdcj) + g^^ dcj)^ + guR dR^ + g^Q dQ'^ , (A3) 

where 

The metric components g^^u in Eq. ()A3|1 are functions of i? and O. 

In deriving Eq. ()A3|1 we have used 

dr dr 09 09 

'"'OROQ^'''OROQ=''^ ^^'^ 

since the coordinate lines of G are orthogonal to the stream foliation surfaces. 

In coordinates {t, R, 9, (p), the Maxwell equation to be solved is again given by Eq. ® , 

but now we have x" = {t,R,Q,(j)) and 

g = 9RR9ee {9ttg4,<t> - gL) ■ (A6) 



Assuming that each magnetic field line lies in one of the stream foliation surfaces, then 

(9?/® OR^ 

Equation ()A7|1 is the generalization of Eqs. (jZj) and (jH])- This assumption should not limit 
our results too much, since the B-component of the magnetic field is not relevant to the 
transportation of angular momentum in the /^-direction. 

Then, for stationary and axisymmetric solutions the Maxwell equation becomes 

^[V^{u'B^-u^B^)]=0. (A8) 



OR 

Making use of UaB"- = and UaU"" = —1, we obtain the solutions to Eq. 

5^ = i-CoUt + ^u^) , (A9) 

= [-Coutu^ + (1 + u^u'^) ^] , (AlO) 

B' = [- (1 + utu') Co + u^u'^] , (All) 
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where Cq = Co(e), ^ = ^(6). From Eqs. and (|XTn|) we have 



B 



R 



u't> 



(A12) 



when \1/ = 0. 

Substituting the solutions ()A9|) - ()A11|) into Eq. ()20|) . we obtain the cj) — R and t — R 
components of the electromagnetic stress-energy tensor 



EM,</ 



R 



EM,t 



(A13) 



where 



(AM) 



The corresponding equations for mass conservation, angular momentum conservation, and 
energy conservation are respectively 



d 



dR i^P-^"") 







_d_ 

OR 
d 



— [v^(p^n,M« + TEM/)] = 



(A15) 
(A16) 
(A17) 



Apparently, when = (i.e., fl^i, = 0), we have T-^ut^ = 0- Then, from Eqs. ()A15|1 and 
()A17jl we have 



dut 
dR 



0, 



(A18) 



i.e., the specific energy E = —Ut is constant on each stream foliation surface. Therefore, 
solutions with constant specific energy exist for the general model considered here. 

Assuming that Ut ~ — m* ~ —1, |m^/(-Ru'^)| ^ 1, and ^ oo a,s R ^ oo, then we have 



R\ 
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\R^oo 



00 



47rpn 



oc 



-g\u' 



00 , (A19) 



R^oo 



if \E' 7^ and {R^ / \/—g) r-^oo is finite or infinitely large. (For references, \/—g oc R^ for 
spherical coordinates, oc R for cylindrical coordinates.) To prevent this to happen, as for 
the model in the text we must have = 0. 
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APPENDIX B: STATIONARY AND AXISYMMETRIC NEWTONIAN ACCRE- 
TION DISKS DRIVEN BY MAGNETIC FIELDS 

The dynamics of a Newtonian, perfectly conducting, and magnetized accretion disk is 
governed by the following equations here we neglect viscosity and resistivity) 



^ + V ■ (p^v) = , (Bl) 

Pm (^^ + V ■ V = -p^ V$ - V (^p + + • VB , (B2) 

— = V X (v X B) , (B3) 

V ■ B = , (B4) 



where pm is the mass density, v is the velocity, $ is the gravitational potential of the central 
compact object, p is the gas pressure, and B is the magnetic field. Equation ()Bljl - ()B4jl are 
respectively the continuity, momentum, induction, and divergence-free equations. 

Let us consider a stationary and axisymmetric Newtonian accretion disk. Similar to the 
relativistic case studied in the text, we adopt spherical coordinates and assume that in a 
small neighborhood of the equatorial plane the velocity and magnetic field are given by 

V = Vrer + v^e^ , B = Brer + B^e^ , (B5) 

where ej {i = r, 6, 0) are the unit coordinate vectors, Vi and Bi depend only on r. Then, 
Eqs. ()BHl - ()B4jl become first-order ordinary differential equations with radius r as the only 
variable. They can be easily solved. 

The solutions to Eqs. (jHH), (jB3l), and (jBH) are 

r^p^Vr = ^m/47r , (B6) 
r {vrB^ - v^Br) = Ci , (B7) 
r^Br = ^M/47r , (B8) 

respectively, where \E'm, Ci, and \&m are integral constants. 

Assume that remains finite as r ^ oo. Then, from Eqs. ()B7|1 and ()B8|) we have B^{r — »• 
oo) = Ci/rvr if Ci 7^ 0. If Vr{r oo) is also finite, then we find that {B'^/inp^)^^^ = 
Cf/'^niVr, where B"^ = B^ + 5?. For the case of accretion disks, it is reasonable to require 
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that Vr = and B'^/Airp^ is finite as r oo, therefore we must have Ci = 0, i.e., 

by Eq. ()B7|) . Equation ()B9|) means that the magnetic field is parallel to the velocity field 
everywhere. 

Substitute Eq. ()B5|) into Eq. ()B2|) . we obtain 



r 2 c/r (ir pm (ir Snp^^r'^ dr ' 

in the r and 0-directions near the equatorial plane, respectively. 
From Eqs. (jB6l), (jHH), and (IbTT|i . we obtain 

1 

Pmt^.^<^ = ^5r5^ + ^ , (B12) 

where \EfL is the constant of angular momentum flux in the radial direction. Substitute 



Eq. (1221) into Eq. ()RT2jl . we obtain 



o2_.„_,2A ^l1 



where L = rv^f, is the specific angular momentum of fluid particles. As r — > cxd we have 
= + i?| = 47rpm'^^5 if -^(^ 00) is unbounded. That is, at large radii the magnetic 

energy density is equipartitioned with the kinetic energy density. Note that this conclusion 

crucially depends on the validity of Eq. ()B9|) . 
From Eqs. ()B9|) and ()B11|) we have 

Substitute Eq. dnTij) into Eq. (inTO|) . we obtain 

^(l^ + t)=-±*. (B15) 

dr \2 J pm dr 

where v"^ = + v'^. For a barotropic equation of state p = p(pm), the solution is 

+ * + = (B16) 

where \E'e is the constant of energy flux, /i(pm) = /q''™ Pnidp is the specific enthalpy. 
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When p = 0, we have h = 0, Eq. ()B16|) then imphes that the specific energy of disk 
particles, (l/2)t>^ + $, is constant. This corresponds to the solutions of constant specific 
energy found in the text. 

From Eqs. ()B6|) . ()B8|) . ()B9|) . ()B12|) . and ()B16|) we can derive the radial momentum 
equation for the case of p = 



(B17) 



.1 - (vl>yi67r2vl/^) (rV)"\ 

Eq. ()B17|1 corresponds to the Newtonian limit of the general relativistic radial momentum 
equation (with \1' = 0). (Note that — ^l/^hi corresponds to /l, — ^|[/167r^^m corre- 
sponds to Cq, — \E'E/^E'm corresponds to /e + 1 and satisfies |\E'E/^E'm| ^ 1-) A differential 
form of the radial momentum equation can be derived from the differentiation of Eq. ()B17|) 
with respect to Vr and r, as we did for the general relativistic case. 

The asymptotic solution for Vr at infinity can be obtained from Eq. ()B17j) . Suppose 
Vr{r —>■ oo) = 0, then Eq. ()B17jl leads to 



1 



2 



\[^E 



Ll-(^^/167r2^^) (rV)" 
as r oo. Since the left-hand side cannot be negative, asymptotic solutions at infinity exist 
only if \l'E/^E'm > 0. Note that — \E'E/^m is just the specific gravitational binding energy 
of fluid particles when p = [see Eq. ()B16|) ]. In this paper we do not discuss a disk with 
negative binding energy (i.e. \l'E/^E'm > 0), thus we will assume \l'E/^l'm = as we discuss 
asymptotic solutions. Then, from Eq. ()B17|) . we have 

*L 1 

- - ^ (B18) 



2 ' 



as r — s> CX3. The corresponding asymptotic solution for is 



.,= H;-2t)"=«f?^£V'\ (B19) 



where we have taken $ = — M/r for the Newtonian potential of a central compact object of 
mass M. The asymptotic solution for pm can be obtained from Eqs. ()B6|) and ()B18|) : 



^ T 



which is a constant \6 



Pm^47r(^) , (B20) 
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The radial magnetic field is always given by Br = \l/M/47rr^ [Eq. ()B8|) ]. The asymp- 
totic toroidal magnetic field can be obtained from Eqs. ()B13|) . ()B19|) . and ()B20|) : i?^ ^ 
±47r (^m/^M) (2M/r)i/2_ Therefore we have 



47r^„ \ ^ 2M 



as r ^ oo. 

All the above asymptotic solutions agree with those listed in Table IJ for the general rela- 
tivistic model. (Note that \l/m corresponds to F^, and \E'm corresponds to 47rCo.) Certainly 
this has to be true since at large distances from the central black hole the dynamics becomes 
Newtonian. 
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TABLE I: Asymptotic solutions corresponding to ^' = and E = 1. Symbols: u'^ , radial com- 
ponent of the four- velocity; Q = u'^/u'', angular velocity; pm, mass-energy density; = BaB°- 
(/9m aiid B'^ are measured by an observer comoving with the flow); L = u^, specific angular mo- 
mentum of fluid particles; E = —ut, specific energy; ratio of electromagnetic energy density 
to kinetic energy density [Eq. (jSH)], measured by a LNRF observer; u"^ is given by Eq. Lh 
by Eq. (|57j) : by Eq. (|60|) {Qq, = 0, /e = —1); ^r, angular velocity of the black hole hori- 
zon. Supersonic/subsonic motion refers to the fact that n^n*" is greater/smaller than c\/ (l — c^); 



UrU'' = c\/ {l- 4) 


defines the fast critical point [Eq. ((51])]. 
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FIG. 1: Solutions for rf — the radius at the fast critical point, corresponding to a magnetized 
flow with constant specific energy accreting onto a Schwarzschild black hole (solid curves). The 
parameter cq, defined by Eq. is proportional to the flux of radial magnetic field. The outer 

boundary of the flow is at the marginally stable circular orbit: tq = r^g = 6rg (dotted line), 
where rg = M is the gravitational radius. The specific energy of fluid particles is equal to that 
of a test particle moving on the marginally stable circular orbit. The solutions for rf consist of 
three branches: A, B, and C, with intersection at the dark point. The dashed curve represents the 
boundary for physical solutions: beyond the dashed curve physical solutions do not exist (see the 
text). (The dashed curve is indeed coincident with the solid curves B and C. To make it visible, 
we have shifted the dashed curve upward a little bit.) 
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FIG. 2: Same as Fig. ^but for vq = lO^rg (dark lines) and tq = oo (gray lines). The specific energy 
of fluid particles is equal to that of a test particle moving on a Keplerian circular orbit at r = vq. 
The vertical dotted line marks the location of the marginally bound circular orbit (rmb = 4rg). 
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FIG. 3: Contour plots of the generation function H66|) versus r and u'^ . The parameters are: cq = 
0.1 Tg, /e = —1, /l = —3. 96522 Tg (corresponding to tq = oo). Then, the fast critical point, which 
corresponds to a saddle point of the generation function, is at (rf,-u[) = (3.98852 rg, —0.053964), as 
indicated by the dark point. Contours are plotted from F = —0.2 to F = 0.2 with step 5F = 0.05. 
The contours of F = are given by the two lines that smoothly pass the fast critical point, the 
one that goes from the right-bottom corner to the left-top corner represents the physical solution 
for the accretion flow. 
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FIG. 4: Upper panel: Solution to the radial momentum equation (|77j) . The parameters are: 
Co = 0.1 Tg, /e = —1, and /l = —3. 96522 rg (corresponding to tq = oo). The dark point shows 
the location of the fast critical point. The dashed line shows the asymptotic solution = — Cg/r^. 
Lower panel: The ratio ^ = [u^ti^ (l — c\) /c^] corresponding to the solution in the upper 
panel. Subsonic motion corresponds to ^ < 1. Supersonic motion corresponds to ^ > 1. The fast 
critical point is at = 1, as indicated by the dark point. 
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FIG. 5: Evolution of the specific angular momentum of the particles in the flow described by the 
solution in Fig. (solid line). The dashed line shows the solution corresponding to a thin Keplerian 
disk joined to a free-fall flow inside the marginally stable circular orbit. The two arrows show the 
locations of the fast critical point (left arrow) corresponding to the solution in Fig. [IJ and the 
marginally stable circular orbit (right arrow). 
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FIG. 6: Evolution of the mass density (solid line) and the ratio of the radial velocity to the 
rotational velocity (dashed line), for the flow described by the solution in Fig. ^ The mass density 
is in units of its value at infinity: p^o = F'^/ {^^-kCq) [see Eq. The dark points show the 

location of the fast critical point. 
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FIG. 7: The ratio of the magnetic energy density to the kinetic energy density [defined by Eq. ()59() ] 
corresponding to the solution in Fig.0I The dark point shows the location of the fast critical point. 
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FIG. 8: The angular velocity of a flow with /e = —1. The parameter cq is alternatively assumed 
to be Irg (solid line) and 10^^ rg (dashed line). The dotted line shows the angular velocity of a 
thin Keplerian disk joined to a free-fall flow inside the marginally stable circular orbit (the dark 
point). The dashed-dotted line represents the asymptotic solution given by Eq. (jHSI)- 
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FIG. 9: Solutions for rf — the radius at the fast critical point corresponding to a flow with constant 
specific energy accreting onto a Kerr black hole (solid curves). Each panel corresponds to a different 
spinning state of the black hole, as indicated by the value of a. The outer boundary of the flow 
is at the marginally stable circular orbit: tq — ^'ms 

(a), as indicated by the vertical dotted line. 
The specific energy of fluid particles is assumed to be equal to that of a test particle moving on 
a Keplerian circular orbit at ro = 

(a). The dashed curve represents the boundary for physical 
solutions: beyond the dashed curve physical solutions do not exist (see the text). 
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FIG. 10: Same as Fig. IHlbut ro = oo. The specific energy of fluid particles is assumed to be equal 
to 1. The vertical dotted line marks the location of the marginally bound circular orbit, r^^^a). 
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FIG. 11: Upper panel: Evolution of the specific angular momentum of a flow with constant specific 
energy accreting onto a Kerr black hole with a/M = 0.95. The dashed curve shows the solution 
corresponding to a thin Keplerian disk joined to a free-fall flow inside the marginally stable circular 
orbit. The two arrows show the locations of the fast critical point (left) and the marginally stable 
circular orbit (right). Lower panel: Ratio of the electromagnetic energy density to the kinetic 
energy density (solid line) corresponding to the solution in the upper panel, which consists of the 
contributions of magnetic field (dashed line) and electric field (dotted line). The dark point shows 
the location of the fast critical point. The circles show location of the horizon. (The unit of the 
horizontal axis is chosen to make the structure near the horizon clear.) 
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FIG. 12: Same as Fig. HU but for 


a/M = 


-0.9. 
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FIG. 13: The specific angular momentum of fluid particles as they reach the horizon of the black 
hole, as a function of cq. Different curves correspond to different spinning states of the black hole 
as labeled. The outer boundary of the flow is at ro = oo so /e = — 1. (At large radii the specific 
angular momentum is cx r^l"^ ^ so L = oo at ro = oo.) 
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